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The outstanding progress in nanostructure fabrication and cooling technologies allows what was 
unthinkable a few decades ago: bringing single-mode mechanical vibrations to the quantum regime. 

The coupling between photon and phonon excitations is a natural source of nonclassical states of 
light and mechanical vibrations, and its study within the field of cavity optomechanics is developing 
lightning-fast. Photonic crystal cavities are highly integrable architectures that have demonstrated 
the strongest optomechanical coupling to date, and should therefore play a central role for such 
hybrid quantum state engineering. In this context, we propose a realistic heralding protocol for 
the on-chip preparation of remotely entangled mechanical states, relying on the state-of-the-art 
optomechanical parameters of a silicon-based nanobeam structure. Pulsed sideband excitation of a 
Stokes process, combined with single photon detection, allows writing a delocalised mechanical Bell 
state in the system, signatures of which can then be read out in the optical field. A measure of 
entanglement in this protocol is provided by the visibility of a characteristic quantum interference 
pattern in the emitted light. 

PACS numbers: 42.50.Wk, 03.67.Bg, 42.50.Dv, 42.70.Qs 


Introduction. Cavity optomechanics m remark¬ 
ably allows to transpose the cavity quantum electrody¬ 
namics features to vibrational quanta, through radiation 
pressure. This field of research has demonstrated a spec¬ 
tacular theoretical and experimental development in the 
past decade and is now envisioned as one of the most 
promising routes to the production of nonclassical states 
of a macroscopic degrees of freedom. Recently a number 
of theoretical investigations have predicted such nonclas¬ 
sical mechanical states to occur in the form of squeezed 
states m\, Fock states ISIS] and photon-phonon ns- 
El or purely mechanical entangled states [T5H21]- Such 
achievements would not only shine light on the founda¬ 
tions of quantum mechanics and the decoherence pro¬ 
cess, but could also provide a long-lived information stor¬ 
age or processing platform useful for potential quan¬ 
tum repeaters [25]. Outstanding experimental results 
have already been obtained in this direction by Lee and 
coworkers [26] who managed, based on the Duan-Lukin- 
Cirac-Zoller protocol to entangle two optical phonon 
modes stored in remote diamond crystals at room tem¬ 
perature. Meanwhile, current cooling technology pSllST] 
is paving the way to nonclassical states stored in acous¬ 
tic phonon modes, which can be pictured as the collective 
motion of a macroscopic number of atoms in the solid. 

The recent prowess in the optimisation of silicon 
[291 [32], indium phosphide [33l [34] or more recently 
diamond-based [35l [36] photonic crystal structures al¬ 
low to combine high quality factors of the cavity modes 
with confined mechanical modes, resulting in record val¬ 
ues of the optomechanical coupling reaching the MHz 
range [29]. Thanks to their small footprint, solid state 
structures lie among the best candidates for the future 
of integrated quantum logic elements that could exploit 
optical and mechanical quanta to code information. In 
this framework, it is desirable to have access to on-chip 


operations that involve single quanta of excitation and 
remotely entangled states [37] of high purity. 

In this article we predict the strong potential of op¬ 
timised silicon nanobeam structures for on-chip entan¬ 
glement of distant mechanical vibrations. We first de¬ 
sign and characterise a promising candidate system based 
on a state-of-the-art photonic-phononic crystal structure. 
Then, we demonstrate the possibility to engineer re¬ 
motely entangled states of mechanical quanta on chip, 
based on a high efficiency heralding protocol [24]. The 
quantum dynamics is described by means of a master 
equation treatment that includes the impact of environ¬ 
mental interactions and finite temperature. 

The structure. We propose the system shown in 
Fig0 which is based on a pair of silicon nanobeams 
of width 529nm and thickness 220nm, at a distance of 
I40nm from each other. The first beam is patterned in 
order to host two defect cavities based on the optimised 
design proposed in Refs. [29l [38] . separated by a distance 
of d = 6.3 jam. The lattice constant of the photonic crys¬ 
tal pattern (cf. Ref. [38] for details) is uq = 436nm. The 
second beam hosts a clogged waveguide section that sup¬ 
ports Fabry-Perot-like modes and serves to couple the 
two cavities in an implementation similar to the one of 
Ref. [39] . Such a structure could be used to remotely cou¬ 
ple cavities that are separated by a distance of hundreds 
of wavelengths. 

To demonstrate the waveguide-mediated optical cou¬ 
pling we model the structure using a commercial-grade 
simulator based on the finite-difference time-domain 
(FDTD) method [40] . By imposing odd, even, or no par¬ 
ity with respect to the x = 0 plane, we can selectively ex¬ 
cite the anti-symmetric (Figj^a)), symmetric (Figj^b)), 
or both modes of the coupled system. The corresponding 
simulated spectra are presented in Figj^c). The modes 
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FIG. 1. 3D model of the system under study, showing the nanobeam photonic crystal (upper) and the clogged waveguide 
(lower), together with a sketch of the excitation and detection conditions. The cavity 1 is driven under pulsed excitation and 
its output is collected by a single photon detector. The yellow arrow highlights the waveguide-induced optical coupling of 
strength Jc between the 2 cavities. The deformation of the cavities illustrate the slight transverse vibrations induced by the 
optical field that are computed in Figs|^d) and (e). 


resonate around the frequency ujc = 27rx 193.45 THz. To 
accurately determine the frequency splitting, we record 
the electric field in the center of one cavity when both 
modes are excited, and observe a clear beating pattern, 
as shown in the inset. From the time interval between 
the two antinodes we compute 2Jc = 27r x 16.9 GHz, 
which is consistent with the spectrum in Figj^c). In 
the absence of the clogged waveguide instead, the simu¬ 
lation results in a splitting of about 4 GHz that would 
dramatically drop for larger cavity separation, thus high¬ 
lighting the role of the waveguide in enhancing the cou¬ 
pling between the two cavities. The quality factors of the 
two modes were computed from the time decay of the 
field after selective excitation and are Q- = 1.14 x 10^ 
and (5+ = 8.0 x 10^, corresponding to linewidths of 
K- = 271 X 169 MHz and = 271 x 243 MHz for the 
antisymmetric and symmetric modes, respectively. The 
linewidths are slightly unbalanced due to the different 
overlap of the two confined modes with the continuum 
of radiation modes outside the beams. Even though 
such high Q-s are experimentally achievable [38], be¬ 
low we shall conservatively assume twice larger values 
of = 27r X 338 MHz and = 27r x 486 MHz in 
view of the extra losses that would be introduced by the 
measurement though input and output channels such as 
tapered fibers m- 

The vibrational eigenmodes of the structure, asso¬ 
ciated to the two defect cavities, were computed via 
finite element modeling m- The localised transverse 
breathing modes, as described in Ref. [29], are found at 
a frequency = 27r x 5.08 GHz. As for the optical 
modes, symmetric and anti-symmetric normal modes 
of the vibrations at the two cavities arise, which are 
characterised by a negligible frequency separation of 
2Jm = 27rx8.8 kHz. The displacement fields of these 
modes are depicted in Figj^d),(e). The associated 
distortion of the cavity is oriented mostly along the 
^-direction, producing a record optomechanical coupling 
of ^ = 27 r X 0.86 MHz [38]. We assume the mechanical 
quality factor Qm = 0.8 x 10^ measured in Ref. m, 



FIG. 2. (a) Anti-symmetric and (b) symmetric normal modes 
of the electromagnetic field in the coupled cavity system. The 
colormap expresses the Re(Ey) component of the electric field, 
(c) Simulated optical spectrum after individual or simultane¬ 
ous excitation of the optical modes. The inset shows the time 
dependence of the electric field, as recorded at the center of 
one cavity when both modes are excited, (d) Anti-symmetric 
and (e) symmetric confined mechanical modes. The colormap 
expresses the amplitude of the displacement field. 


where the structure was coupled to an acoustic radiation 
shield. As we will see, the readout of the mechanical 
entangled state requires that the vibrational modes 
are nondegenerate, in order to introduce a relative 
phase factor between the two optomechanical coupling 
processes, that reveals the quantum superposition of the 
two mechanical modes through an interference pattern 
in the optical readout m- Different frequencies for the 
mechanical modes of the two cavities will naturally arise 
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due to imperfections in the nanofabrication process. 
We have therefore studied the impact of the structural 
disorder on the resonant mechanical frequency of one 
cavity, by introducing a random, Gauss distributed 
deviation of the holes dimensions and positions with a 
standard deviation of 1 nm. Over 100 simulated disorder 
realisations we found a typical 1% fluctuation of 
that shall therefore be assumed from now on namely 
^m,2 = 1.01 X 0^,2 = 27r X 5.13 GHz. 


Quantum model. The system can be schematically 
summarised as in FigJ^a) namely as two coupled, single¬ 
mode cavities (di^2 operators), each of them optomechan- 
ically coupled to distinct mechanical modes (61,2 opera¬ 
tors). The system Hamiltonian reads 


i~L — ^ ^ 

i=i,2 

-h Jc (^b\d2 + + F (t) -h 


( 1 ) 


We introduce an optical write and readout protocol anal¬ 
ogous to the one we recently adopted in the case of a sin¬ 
gle optical cavity m- The system is driven first by a clas¬ 
sical source of amplitude F{t) and frequency ujl on the 
cavity 1 that stands for our input-output mode (see W 
in Figj^c)). It is instructive to rewrite on the basis of 
the symmetric and antisymmetric modes following a± = 
(di ± 0 . 2 )IV 2 and b± = (bi ± 1 ) 2 )I which recasts the 
optomechanical contribution iig = g^j dtdj(6t + hj) as 

^ + ald_^ (bl_ + 6+^ (2) 

+ -^= . 

It describes two cavity modes of resonance uj± = ujc^ Jc 
(i) optomechanically coupled to the same vibrational 
mode 6+ of frequency 11+ = (Hi -1- H2)/2 — Hi^2 (first 
term) and (ii) coupied together through the absorbtion or 
emission of a phonon of the mode (second term). The 
iatter process, that wouid aiiow an exchange of quanta 
between the d+ and d_ modes, is off-resonant in the con¬ 
ditions we consider where co’+ — uj- 3Hi^2 and can be 
discarded. The cavity spectrum shown in Figj^c) con¬ 
sists of two main resonances at cjc di J each surrounded 
by sidebands formed by the Raman interaction with the 
5+ mode (dashed fines). Note that (ii) makes the d+ 
sidebands visibie in the d_ spectrum and vice versa. 

If the cavity 1 is driven on the resoived d+ antistokes 
sideband, i.e. setting ujl = uJcF J F H+ (see W arrow in 
FigJ^c)), the first term of Eq.Q then produces the two 
mode squeezing [4] 


^ ^ (d\_b\_ + , (3) 


The d_ contribution is negiigibie since it is shifted by 
— (2J -1- H+) from the pump frequency which makes it 




FIG. 3. (a) Schematic representation of the model system, (b) 
Write and readout pulses sequences for the interference recon¬ 
struction. (c) The computed power spectrum of the symmet¬ 
ric and antisymmetric cavity modes, revealing the stokes and 
antistokes sidebands. 


strongly nonresonant. The interaction (§ therefore cre¬ 
ates a pair correlation between the symmetric cavity fluc¬ 
tuations and mechanical modes, resulting in a squeezed 
state of the mode 6+. Under continuous wave driving, 
one could therefore achieve a steady state hybrid continu¬ 
ous variable entanglement between the d+ and 6+ modes 
m- However, the goal here is to produce a mechanical 
excitation in a linear superposition of two Fock states, 
each corresponding to one quantum of vibration in one 
cavity. This requires combining pulsed excitation with 
postselection. Indeed, after the write pulse, the detec¬ 
tion of a single photon emitted by any of the two cavities 
heralds the presence of a single phonon [SlEl] stored in 
the mechanical mode 5+, i.e. shared between the remote 
modes bi and 62. 

Assuming zero temperature, the system lies initially 
in its vacuum state \ 2 po) = |0)^^ G |00)^ G) |0)c2 5 where 
the indices Cj and m refer to the cavities and mechanical 
states respectively. The write laser tuned to the Stokes 
process of the d+ mode ideally results in the state 

iig I'lpo) = Coo|0)^^ ® |00)„ (8) |0)^^ 

+ cio|l),,®(|lO)„ + e*‘^|Ol)j0|O+ (4) 

+ coilO),^ ® (|10)„ + e^^^lOl)^) ® |1),^ + ... 

Here (p = (po (H2 — Hi)t where (po is a constant phase 
shift depending on the details of the excitation scheme, 
and the coefficients Cij are determined by the driving 
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field. In particular one expects that, in the limit of low 
driving field, the components of this state corresponding 
to two- and more photons are negligible. A single photon 
detection from the cavity 1 retains only the second term 
of Eq. and the state of the system, after detection, 
reduces to 

\^w) = c'lolO),, ® (|10)„ + e*‘^'|01) J ® |0),^ (5) 

which is nothing but a maximally entangled mechanical 
Bell state \i>H)m = (|10)m + exp(V)|01)„)/v^ when 
the cavity modes are traced out. 

The nonclassical character of the mechanical quantum 
state can only be indirectly assessed through the optical 
field. For this, the information stored in the phonon state 
has to be transferred back to the cavity field. This is 
achieved thanks to the anti-Stokes process induced by 
a second pulse of frequency ujl = cjc + T — onto 
the cavity 1 (see R arrow in FigJ^c)). Such excitation 
condition is resonant with the beam splitter interaction 

m 

~ ^ («+b + a+^+) , (6) 


governed by the following set of equations 
ia± = -]- ^a/ 2 Re {I3±) — 

+ ^V^Re (^zp) Qfip H (8) 

ip± = 

+ ^ (9) 

that are written here in the frame rotating at the cavity 
frequency uJc and on the basis of symmetric and antisym¬ 
metric modes a± and P± to account for their unbalanced 
losses i<i±. We however consider equal mechanical losses 
7 ± = 27r X 3.75 kHz which is valid given the large time 
scales involved. The dynamics of quantum fluctuations 
is treated via a master equation for the density matrix 
in the displaced reference frame set by the classical held 

m 


dp 

dt 


—i 




i=± 


( 10 ) 


\ F + 'ljnth'£>[Sb]]p 

i=± 


that swaps the 6+ phonon with an d+ photon. When ^ 
is applied to the written state one obtains 

IV’ij) = 4 'o |07 0 (l + |11)„ ® 107 

+ c'io|iy ® (l + e'‘"")|00)n^®|0)c. (7) 

+ <^0110)01 0 (1 + ) | 00 )^ 0 |l)g 2 + ... 

The corresponding cavity 1 intensity Xi = 
(V’iil api IV’ii) = kioRl + cos ((/?")], where 
p'' = Po depends on the delay from the 

write procedure and presents a full contrast interference 
pattern providing means to reveal the entangled nature 
of the mechanical state. Note that one could equally 
detect the cavity 2 held and obtain the same result. 
The expected interferences can be reconstructed by 
integrating the emission received by a single photon 
detector over several runs where the delay between the 
write and readout pulse is gradually increased as shown 
in FigJ^b) and detailed below. 

The simple argument based on pure states that was 
presented above can be supported by a numerical study 
of the full system dynamics, that we carried out using the 
theoretical tools to model open quantum systems [43] . 
In this analysis, in order to be able to account for the 
sizeable classical held component induced by the driving 
laser, we separated the held into the sum of its classical 
and quantum-fluctuation components according to d± = 
a± + Sd± and b± = /3± Sb±. The classical evolution is 


where nth is the mean thermal phonon occupation and 
Vid] p = dp + pd^d — 2dpd^ describe the dissipations in 
the environment for each mode. 

Results and discussion. The write procedure is con¬ 
ducted with a Gaussian excitation at a time tw = 30 ns, 
of amplitude Aw = 10^K-, standard deviation at = 3.85 
ns and frequency ujl = uJc-\- J on the cavity 1. Note 
that one could equally excite the antisymmetric mode at 
= cJc — J + as prescribed by <§• The driving pulse 

results in a large classical held amplitude for the cavity 
mode. This acts as a source for the mechanical held 
and produces a sizeable phonon occupation in both me¬ 
chanical modes, thanks to the optical coupling between 
the two cavities. In this way, the normal modes d+ and 
6+ are excited without needing to drive the held in the 
second cavity. Figures [^a) and (b) show respectively 
the classical and quantum fluctuation occupations. We 
observe that the classical components vanish along with 
the input pulse. Starting from a time tq 25 ns (vertical 
black line), quantum fluctuations dominate over the clas¬ 
sical components and decay according to the respective 
lifetimes. The quantity tq deflnes a time horizon past 
which the system evolves freely and can display quan¬ 
tum correlations. Importantly, a compromise must be 
made on the pulse duration, which should be set (i) long 
enough to accurately excite the sideband without leak¬ 
ing too much on the main cavity frequency and (ii) short 
enough to allow for a sufflciently large occupation of the 
quantum fluctuations at time tq. 
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FIG. 4. Write procedure: Log scale (a) classical and (b) fluc¬ 
tuation occupations versus time. Blue and dashed-red lines: 
cavity mode 1 and 2; Yellow and dashed purple lines: mechan¬ 
ical mode 1 and 2. The vertical black line marks the time tq, 
starting from which quantum fluctuations become dominant, 
(c) The computed heralding rate (t) versus the concurrence 
C{t). (d) Absolute value of the heralded mechanical density 
matrix at the Cmax time (vertical red dashed line in (a) 
and (b)) conditioned on the detection of a single photon out 
of cavity 1. 


The heralding at a time tpi > tq is modeled via the 
projection of the density matrix onto the single photon 
subspace [Sin of cavity 1 


.( 1 ) _ Vip{t^)Vi 

Pn 


-1 


Tr 


Vip{tu) 


A = |l)e,(l|c,^£ 


> Im2 ^ ^C2 


( 11 ) 

( 12 ) 


We consider in Eq. 0 the full density matrix of the sys¬ 
tem, obtained by including the classical contributions as 
coherent states \a±) and \/3±) and applying a multimode 
displacement operator 'D{a±,/3±) to p. The heralding 
efficiency - therefore the degree of mechanical entangle¬ 
ment - is computed through the concurrence C [45] of 
the reduced density matrix pm = TTg,^^g,^{ppi) restrained 
to the 0 and 1 phonon subspaces. It displays a large 
maximum value of Cmax = 0.95 from = 30 ns (see 
vertical dashed lines in FigQa),(b)). Corresponding to 
this value, a mechanical Bell state storing delocalised 
phonon is formed, as revealed by the visual shown 
in FigQd). The imbalance in the diagonal elements is 
a consequence of having driven a single cavity only. We 
have computed a fidelity = 0.97 to the (|10) + |01))/v^ 
Bell state. Obviously, the heralding time would be ran¬ 
dom in an experiment and the detector should target the 


t > tq time window where C ^ Cmax, while avoiding to 
monitor earlier times, for the procedure to succeed. The 
entangled state lives as long as the mechanical lifetime 
Tm = 27r/7i^2- Given the very high quality factor charac¬ 
terising nanobeam cavities, combined to the GHz-range 
of the frequency of the mechanical modes, one could ex¬ 
perimentally reach = 277 ms [38] which is promising 
for quantum repeater technology [25] . When aiming at a 
concurrence C > 0.9, the largest heralding rate achieved 
in our simulations was upi = - 3 X 10“^ MHz, 

where k,i ^2 = + ^-)/2, as shown in the figure j^c). 

Finally we point out that, in addition to a large value of 
the simulated concurrence, in order to assess our target 
entangled state it is required that the two- and higher- 
photon terms in the full mechanical density matrix be 
negligible as it is the case here. This is simply achieved 
by controlling the amplitude of the driving field. 



t-At|^(ns) T(K) 


FIG. 5. Readout procedure: Logarithmic scale (a) classical 
and (b) fluctuation contributions to the state occupations ver¬ 
sus time, as computed for a given realisation of the readout 
procedure. The inset of panel (b) shows the mechanical oc¬ 
cupation in linear scale, (c) The held intensity in cavity 1, 
averaged over several runs, as a function of the readout tn 
and shifted evolution time t — Atn where Atn — t^^ — t^^ for 
the run number j. The inset shows a (normalised) cut along 
the vertical red dashed line, from which the visibility V can be 
inferred. The horizontal white dashed line corresponds to the 
(a) and (b) evolutions where the intensity is at a maximum, 
(d) Simulated temperature dependence of the negativity and 
the fringe visibility. 


The readout step involves a second input pulse hav¬ 
ing the same width as the write pulse, and amplitude 
J[r = 5 X 10 ^delayed by AT = tji — tw The cen¬ 
tral frequency of the pulse is tuned to the anti-Stokes d+ 
sideband ujl = cjc + T — and aims at transferring 
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the information about the phonon mode back to the op¬ 
tical field. As discussed previously, the nondegeneracy of 
the mechanical modes AQ = Q 2 — MHz, result¬ 

ing from the sample imperfections, introduces a relative 
phase factor (p = (po AClt between the modes which 
is associated with an interference pattern in the emit¬ 
ted light. The latter is revealed by varying the write 
and readout pulse delay AT and integrating the corre¬ 
sponding signals (see FigJ^b)). Figures [^a),(b) show the 
classical and fluctuations contributions, respectively, to 
the occupations in a realisation where the cavity intensity 
reaches a maximum (dashed white line in panel(c)). Once 
again, quantum fluctuations dominate after the driving 
pulse has decayed (vertical black line). Starting from 
this time, the quantum interference should be recorded. 
From the panel (b) inset one can see the 6 + mode being 
emptied by the optical field as revealed in the 61 and 62 
simultaneous population decay. 

Varying gradually the delay AT between the write and 
readout pulses and specifically Ir reveals the predicted 
interference pattern of period T = 27r/Af^^ 125 ns, 

as shown in FigJ^c). The field in cavity 2 obviously dis¬ 
plays a similar behavior (not shown). Note that to ob¬ 
tain a vertical pattern in the plot (instead of an oblique 
one), the readout pulses have been shifted so that the x- 
axis shows the the shifted evolution time t — AIr where 
= t^R —t^R j is the run index. This corresponds 
to what would be obtained varying the phase po of the 
write pulse instead of AT. A normalised cut along the 
vertical axis (dashed red line), taken at fixed delay from 
tR, is plotted in the inset. It shows a fringe visibility as 
large as V ^ 0.92. The interference carries a signature 
of the relative phase between the two confined phonon 
modes. This can be related to the amount of phonon 
entanglement. In particular m, a fully separable state 
would result in a visibility Vc = 0.5. This value sets a 
lower bound, below which the entanglement may still be 
present in the system, but cannot be inferred from the 
interference pattern. The average cavity occupation at 
this value of Ir amounts to ni ^ 5 x 10 “^, which cor¬ 
responds to a photoemission rate of = 0.2 MHz. This 
value is a consequence of the large driving that temporar¬ 
ily displaces the system to a coherent state [46] until its 
quantum nature is restored. 

At zero temperature we obtained a Bell state with 
high fidelity, which could be characterised by the con¬ 
currence, if restricting the analysis to the subspace with 
0- and I-phonon states. However, at finite temperature 
a non-negligible thermal occupation may result in a siz¬ 
able contribution from states with n> 1 phonons. In this 
case, an entanglement measure appropriate to the larger 
Hilbert space is needed. The concurrence was recently 
extended to 4-dimensional spaces m- Its computation 
however is a demanding numerical task for n > 4 |48l [49] 
and is therefore not appropriate in the general case we 
want to consider here. The Bell states naturally suggest 


to deploy an entanglement monotone derived from posi¬ 
tive partial transpose criterion, namely the negativity ISO] 
V= 11/3^^'^ 111 — I, where is the partial transpose of 
prn and llpm’^lli its trace norm. The above argument 
using quantum states can be extended to mixed states, 
in order to estimate a lower bound on the visibility re¬ 
sulting from an entangled state, assuming in particular 
that the mechanical modes initially lie in the product 
of 2 thermal states. In the supplementary information 
[51] we analytically derive the Heralded mechanical den¬ 
sity matrix at finite temperature, and recover the lower 
bound of V ^ 0.5 to the fringe visibility produced by an 
entangled state. 

Besides, the impact of state-of-the-art cooling tem¬ 
perature T m on our protocol was computed by 
direct master equation simulations, and the results 
are presented in FigJ^d), where Af and V versus T is 
plotted. We see, as expected, a clear correlation between 
JV and V. They both start decreasing from T 0.05K. 
Entanglement survives in the temperature range we 
considered, and the visibility reaches the lower bound Vc 
at Tc ^ 0 . 2 K. Below this temperature, our protocol can’t 
assess the presence of entanglement, and a full quantum 
tomography of the state would be instead required. 
Finally the impact of the pure dephasing imposed by 
other mechanical resonances was analysed in Ref.[24| 
and was shown to be negligible in the photonic crystal 
parameter range. 

Conclusion. We have proposed a realistic photonic 
crystal architecture based on coupled nanobeams and 
computed its optomechanical properties. We have 
shown that this structure should allow to achieve 
on-chip remote entanglement between the localised 
acoustic phonon modes, mediated by the optical field, 
with no need for an additional beam splitter as in typical 
entanglement protocols [ini US]. The protocol that we 
propose achieves high heralding rates and is still effective 
at typical cooling temperature for Silicon-based systems. 
The characteristic long mechanical lifetime makes 
the present system very promising for demonstrating 
macroscopic entanglement and for the future of hybrid 
quantum logic circuits. 
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